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Thermodynamic data on the n-paraffins from n-C 6 Hi 4 through n-Ci8H 38 have been 
used to obtain values for the specific heat, entropy, enthalpy, and Gibbs free energy of a 
large, ideal CH 2 -chain crystal from to 420 °K and of the liquid above 200 °K. Analytical 
expressions are given for the properties of the crystal and liquid above 200 °K. For the 
crystal, a modified Einstein function was used to adjust the melting temperature to 414.3 °K. 
Values between 975 and 1025 cal/mole for the heat of fusion were found to be the ones most 
consistent with the data. Comparison of the results with polyethylene data shows reasonable 
agreement. 



1. Introduction 

This paper concerns the empirical determination 
of accurate values for the thermodynamic properties 
of ideal CH 2 -chain crystal 1 and liquid in the temper- 
ature range from to 420 °K and 200 to 420 °K 
respectively, using published data on pure n-paraffins. 
This analysis is intended to circumvent the inherent 
difficulties involved in measuring polyethylene 
directly since (so far) even the most carefully 
crystallized polyethylene samples are not large, well 
defined CH 2 -chain systems, and measurements neces- 
sarily reflect the nonnegligible effects of chain folds 
and "amorphous" inclusions. Wunderlich [l], 2 
realizing this problem, has extrapolated specific 
heats to 100 percent crystallinity to obtain the 
crystalline CH 2 contribution to the specific heat. 
This method is presently limited in accuracy by a 
scarcity of data and variations in the chemical 
structures (e.g., degree of branching and molecular 
weight) which exist between samples of different 
crystallinity. Even ideally, if one had sufficient 
data on several polyethylene samples which differed 
in percent crystallinity, the extrapolated results 
might still reflect the presence of chain folded 
lamellar surfaces, unless the lamellar thicknesses were 
properly taken into account. 

The use of n-paraffin data to determine CH 2 prop- 
erties offers the advantages that the 7i-paraffins can 
be made highly pure and crystallized to a high 



i This means a crystal, large in all dimensions, of very long CH2 chain mole- 
^ cules, which contains no chain ends or branches and only an equilibrium number 
of defects. 
2 Figures in brackets indicate the literature references at the end of this paper. 



degree of perfection. Very carefully measured spe- 
cific heat data are available for pure /^-paraffins from 
u-CqR u through n-Ci 8 H 38 (see following section). It 
has been reported that the molar specific heat (C p ), s 
enthalpy (/Z), entropy ($), and Gibbs free energy 
(F) of the solid n-paraffins appear to be linear with 
chain length at constant temperature and crystal 
phase [2-4]. An obvious interpretation of this 
observation is that each CH 2 unit in the chain can 
be thought of as contributing a certain number of 
vibrational modes to the crystal without changing 
the frequency spectrum, and that this contribution 
is independent of the length of the chain in which 
the unit is located, so long as the local crystal 
structure remains unchanged. Accordingly, the 
properties of a CH 2 unit should be the same in a 
large ideal — CH 2 — crystal as in the ^-paraffins. 
The major limitation to the n-paraffin determination 
of the crystalline CH 2 properties are the deviations 
from linearity with chain length which become appre- 
ciable at temperatures near the n-paraffin melting 
points. 

For the liquid phase, there is a slight but real 
departure from linearity which will necessitate special 
treatment in the determination of the liquid phase 
properties. 

The approach used below will be: (a) To establish 
the best values for C p , H, S, and F of a — CH 2 — 
unit in both the liquid and solid phases over as wide 
a range of temperatures as possible, (b) to establish 
a consistent analytical representation of the above 

3 Actually, the specific heat data refers to C 8a td, the specific heat under the 
saturated vapor pressure of the material. For this work, the difference between 
C p and <7 Ba td for the solid and liquid phases of the rc-paraffins is negligible. 
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data from 200 to 420 °K, and (c) to compare the 
results with actual data on polyethylene and examine 
l he discrepancies between the predicted and observed 
properties. 

2. Source of Data 

The 71-paraffin data used in this study were 
furnished by the United States Bureau of Mines, 
Bartlesville, Okla., in the form of tables of values 
for C p , H-H , S, and F-H (H is the enthalpy at 
°K). Their work on ti-C 6 H 14 , w-C 7 H 16 , and 
n-CgRis through n-C 16 H 34 is presented, together with 
details of sample purity, measurement techniques, 
and thermodynamic calculations, in references [5, 6, 
and 7] respectively, and their work on w-Ci 7 H 36 and 
7?-Ci 8 H 38 is as yet unpublished. The data on 
n-C 7 H 16 , 7i-Ci 7 H 36 and n-Ci g H 38 covered the tempera- 
ture range from to 370 °K, and data for the remain- 
ing paraffins extended from to 320 °K. 

3. CH 2 Structures 

The even n-paraffins 4 from C 6 H 14 through Ci S H 38 
are triclinic at all temperatures up to the melting 
point. The odd n-paraffms 4 from CnH 24 through 
Ci 7 H 36 are orthorhombic at all temperatures up to 
within 15 degrees of their melting points but under- 
go a solid-solid transition prior to melting. The 
orthorhombic CH 2 structure of the odd paraffins 
appears to be identical to the CH 2 structure of poly- 
ethylene, and is assumed to be identical to the struc- 
ture of the ideal CH 2 -chain crystal whose properties 
we are predicting. Although one would not expect 
the data for the triclinic structure of the even n- 
paraffins to be relevant to this work, it was found 
that the data obtained from the even /i-paraffins 
varied only slightly from that of the odd 7i-paraffins. 

In the liquid state one expects and finds no even- 
odd differences among the paraffins. The absolute 
entropies of all 13 paraffins examined here, when 
plotted against chain length, lie on the same, nearly 
linear curve and the enthalpies and free energies 
show a similar behavior if one takes into account 
the differences in zero point enthalpies between 
the different solid structures. For instance, the 
enthalpy which had to be added to F and H for the 
even paraffins in order to have the liquid phase 
values of F and H for both even and odd paraffins 
lie on the same nearly linear curve, when plotted 
versus n, was found to be 370±10 cal/mole of par- 
affin molecules. This zero point energy difference is 
independent of chain length since it results en- 
tirely from the differences in packing of the chain 
end groups, with the triclinic phase allowing the 
lowest energy packing. 5 

Whereas the solid phase of very linear, high molecu- 
lar weight polyethylene contains a significant volume 
of chain folds which cause deviations from the ideal 
crystalline CH 2 specific heats, the liquid state should 
present no such difficulties. Specific heat measure- 
ments should not be able to detect the small per- 

< Even 7i-parafiins and odd n-paraffins are terms used here to mean w-parafFms 
containing an even or odd number of carbon atoms in the molecule. 
6 See reference [3] for a further discussion of the above points. 



centage of chain ends and branches present in the 
liquid. Therefore, the specific heats and absolute 
entropies calculated here from the liquid ^-paraffin 
data should be close to the actual measured values 
for polyethylene, and a comparison of predicted and 
experimental properties of the liquid provides an 
experimental test of this work. The liquid paraffin 
data should not, however, be expected to reflect the 
presence of a glassy state even well below the tem- 
perature where the ideal CH 2 -chain liquid might 
become a glass. 

4. Analysis of the Paraffin Data 

Table 1 contains a sample of the thermodynamic 
data as received for the solid and liquid phases at 
representative temperatures. 

For the solid at a given temperature, the differ- 
ences between the thermodynamic values of consec- 
utive even and consecutive odd compounds was com- 
puted and these differences were averaged over the 
whole range of compounds for which data were 
available at that temperature. The averages for 
the even and odd paraffins were found to be very 
similar, and the combined average for both even and 
odd compounds was divided by two to obtain the 
desired thermodynamic value per mole of CH 2 units. 

At the higher temperatures, where the shorter 
paraffins are liquid, the accuracy of the computations 
decreases due to the decreased range of data and the 
premelting 6 increase in the specific heat that begins 
about 50 degrees below the melting point and which 
can be seen in the data at 200 °K in table 1. It is 
probable that above 200 °K, the values calculated 
by the above procedure will be somewhat smaller 
than the correct values. 

For the liquid phase data, the range of tempera- 
tures amenable to computation was limited at low 
temperatures by the decrease in the number of liquid 
samples and at high temperatures (370 °K) by the 
termination of experimental measurements. The 
most reliable values are those around 300 °K where 
most data exist. 

The differences between the thermodynamic 
values for consecutive even and odd paraffins are 
not quite as constant for the liquid as for the solid 
as can be seen in table 1, and an appropriate ex- 
trapolation of these differences to n= °° was con- 
sidered to be preferable to using a simple average. 
When plotted versus 1/n, the differences for the liquid 
were found to be linear within experimental error. 
For F } H, and S, the best linear extrapolation to 
l/n = was graphically determined, and the resulting 
intercept was chosen as the desired thermodynamic 
value. For C p , this \\n extrapolation gave values 
which exhibited excessive scatter. It was found 
that this scatter could be reduced by integrating 
the G v data, extrapolating the integrated values to 
l/7i = and then differentiating the extrapolated 
integrated values to obtain the desired extrapolated 
C p values. In practice, this process merely involved 



e We are not here concerned with the distinction between homophase premelt- 
ing and impurity induced heterophase premelting. 
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Table 


1. Sample 


data and calculations 
















Number of carbon atoms in molecule, n 


Thermodynamic 
function 


6 


8 


10 


12 


14 


16 


18 


Even 
avg. 


7 


9 


11 


13 


15 


17 


Odd 

avg. 


Net 
avg. 
per 
C Hj 


Graph, 
extrap. 

to 
l/n=0 




Solid phase at 100 °K 


Cp(cal/mole °K) 




24.08 | 28.60 | 33.16 | 37.72 42.26 46.86 
4.52 4.56 4.56 4.54 4.60 

20.47 | 23.75 27.05 | 30.36 33.67 36.98 
3.28 3.30 3.31 3.31 3.31 

1207 1414 1622 1829 I 2036 I 2243 
207 208 207 207 207 

840 961 1083 1207 | 1331 1454 
121 122 124 124 123 








31.17 35.79 40.27 44.84 
4. 62 4. 48 4. 57 

25. 55 1 28. 89 1 32. 16 35. 55 
3. 34 3. 27 3. 39 

1525 1735 1940 ' 2151 
210 205 211 

1031 1154 1276 | 1403 
124 122 127 








C P (n+2)— C p (n) 




4.56 






4.55 


2. 28 




Sfcal mole °K) 










S(n+-2)— S(n) 




3.30 






3.33 


1. 66 




//"(cal/mole) 










H(n+2)—H(n) 




207 






209 


104 




--/'(cal /mole) . 










-F(n+2)+F(n) 




123 






124 


62 
















Solid phase at 200 °K 


Cp(cal/mole °K) 




39.18 45.86 53.01 | 60.29 67.63 74.98 
6.69 s * 7.15=* 7.28 7.34 7.35 

41.77 48.99 | 56.31 | 63.68 71.02 78.45 
7.22* 7.32 7.37 7.34 7.43 








50.67 57.78 64.94 | 72.31 
7.11a 7.16a 7.37 

53.15 60.60 67.90 | 75.35 
7. 45 7. 30 7. 45 

5615 | G432 7234 8048 
817 802 814 




1 


C p (n+2)—C P 0)-_- 




7.32 






7.37 


3.67 . 


S(cal/mole °K) 












S{n+2)— S{n) 




7.36 






7.39 


3. 69 




//(cal/mole) .... 




4364 1 5152 | 5954 | 6764 | 7568 
788" 802 810 804 8 


8385 








H(n+2)—II{n) 




7 


808 






811 


405 


. 


-/<'( cal/mole) _ _ 




3992 1 4646 5308 5973 ! 6637 7304 
654" 662 665 664 667 






5016 5687 6346 
671 669 6 


7022 




-F(n+2)+F(n) 




664 






76 


672 


334 




















Liquid phase al 300 °K 






C p (cal/mole °K)__._ 


46.86 | 60.90 | 75.38 90.09 105.01 | 120.07 
14.04 14.48 14.71 14.92 15.06 

71.05 1 86.67 | 102.27 117.87 1 133.41 1 148.86 
15.62 15.60 15.60 15.54 15.45 

11301 1 14807 1 18328 21854 25368 | 28871 
3506 3521 3526 3514 3503 

10014 11193 12353 13506 14654 1 15788 










82.68 1 97.47 | 112.53 | 127.87 
14.79 15.06 15.34 

110.02 12fi.Jffl 141.11 156 .68 








C P («+2)— C„(») 




14.64 






15.06 


7.4 


8.0 


S(cal/mole °K). 


78.86 


94 50 




S(n+2)—S(,i, 




15.57 


15.64 15.52 15.57 15.52 15.57 

12681 | 16204 1 19719 23240 | 26783 30320 
3523 3515 3521 3543 3537 

10977 1214K 13287 144M7 IfififiO 16683 


15.58 


7.79 


7.74 


//(cal/mol^... 




H(n+2)— H(ri) 

- F(cal/mole) 




3514 


3526 


1760 


1770 


-F(n+2)+F(n) 


11 


79 1160 1153 1148 1134 




1155 


11 


12 1139 1150 1113 1133 


1141 


574 


547 



1 These values are low because of premelting and were omitted from the calculations. 



taking differences between the consecutive experi- 
mental II values in table 2 and dividing by 10°. 
For example, C v (per mole of CH 2 units) at 305 °K 
was taken to equal II (per mole of CH 2 units) at 
310 °K minus U (per mole of CH 2 units) at 300 °K 
divided by 10°. The same procedure was followed 
(though mathematically, less precisely) using the 
TS data, in order to increase the number of available 
C p values. 

The results of the above determinations are listed 
in table 2 and are represented graphically as closed 
circles in figures 1 through 4. The large number of 
digits given in the calculated value columns in table 
2 is not intended to indicate significant figures but 
is included only for internal consistency. 

5. Determination of Analytical Expressions 

Linear least squares fits of the specific heats for the 
solid above 150 °K and specific heats for the liquid 
above 200 °K as a function of temperature were 
carried out. The calculated standard deviations 
were 0.15 c l/mole-°K for the solid (39 data points) 
and 0.84 cal mole-°K for the liquid (36 data points). 
Integrations of C v and C p /T were then carried out 
and the average values of the constants of integra- 



tion were determined from the entropy and enthalpy 
data. The resulting functions fit the II and S data 
for the solid above 200 °K with a standard deviation 
of 2 cal/mole and 0.01 cal/mole-°K, respectively, 
and for the liquid above 200 °K with a standard 
deviation of 6 cal/mole and 0.02 cal/mole-°K, 
respectively. The Gibbs free energy was computed 
directly from H and S and was found to fit the free 
energy data above 200 °K with a standard deviation 
of 4 cal/mole for the solid and 9 cal/mole for the 
liquid. 

By using the calculated analytical expressions to 
extrapolate up to the temperature where the free 
energies for the solid and liquid become equal, it was 
found that the predicted melting point was 397.4 
°K and the predicted heat of fusion was 1132 cal/ 
mole. Since the predicted melting point is 4 percent 
below its most probable value of 414.3 °K [4] and 
since the predicted heat of fusion seems unreasonably 
high, the assumption of a linear specific heat for the 
solid above 200 °K was concluded to be inadequate 
for the representation of the known data. Hence, 
a correction to the specific heat for the solid was 
sought that would decrease the free energy of the 
solid by a sufficient amount to increase the melting 
point to its previously determined most probable 
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Table 2. Values of the molar specific heats (C p ) enthalpies 
(H — H ), entropies (S) and Gibbs free energies (F — H ) of 
ideal CH 2 chain crystal and liquid as derived from n-paraffin 
data (experimental values) and as calculated using the equa- 
tions in table 3 





C p (cal/mole 


H-H (cal/mole) 


S(cal/mole 


-(F 


-Ho) 


Tem- 





K) 








3 K) 


(cal/mole) 


pera- 


















ture 




















exp. 


calc. 


exp. 


calc. 


exp. 


calc. 


exp. 


calc. 




Solid phase 


°K 


















10 


0.025 
.15 
.39 
.70 
1.03 
1.34 
1.63 
1.88 
2.10 
2.28 
2.45 
2.60 
2.75 
2.89 
3.05 
3.16 
3.28 
3.40 
3.52 
3.67 


3.02 
3.17 
3.32 

3.48 
3.65 
3. 8572 


0.06 

.85 

3.50 

9.00 

17.8 

29.5 

44.6 

62.0 

82.0 

104 

128 

153 

180 

208 

237 

268 

300 

334 

369 

405 




0.01 

.05 

.16 

.32 

.51 

.73 

.95 

1.19 

1.43 

1.66 

1.88 

2.10 

2.32 

2.52 

2.73 

2.92 

3.12 

3.31 

3.51 

3.69 




0.02 
.30 
1.34 
3.7 
7.8 
14.0 
22.4 
33.2 
46.2 
61.7 
79 
99 
121 
146 
171 
200 
230 
263 
297 
334 




20 








30 








40 








50 








60 








70 








80 








90 








100 








110 








120 








130 








140 








150 








160 








170 








180 








190 








200 


405. 88 


3. 7201 


338. 14 


210 


3.81 


4. 0378 


441 


445. 35 


3.87 


3. 9127 


371 


376. 30 


220 


3.95 


4. 2224 


479 


486. 64 


4.05 


4. 1046 


410 


416.39 


230 


4.05 


4. 4107 


519 


529. 78 


4.23 


4. 2963 


453 


458. 39 


240 


4.16 


4. 6028 


562 


574. 85 


4.41 


4. 4881 


496 


502. 31 


250 


4.32 


4. 7975 


603 


621. 84 


4.58 


4. 6799 


541 


548. 15 


260 




4. 9936 


647 


670. 76 


4.77 


4. 8717 


590 


595. 90 


270 




5. 1924 


691 


721. 70 


4.94 


5. 0640 


643 


645. 59 


280 




5. 3909 


742 


774. 57 


5.13 


5. 2562 


695 


697. 18 


290 




5. 5838 


790 


829. 37 


5.30 


5. 4484 


748 


750. 67 


300 




5. 7898 
5. 9909 
6. 1851 
6. 3836 

6. 5795 
6. 7734 




886. 34 
945. 19 
1005. 91 
1068. 84 
1133.64 
1200. 32 




5. 6416 

5. 8345 

6. 0272 
6. 2208 
6. 4142 
6. 6075 




806 15 


310 






863. 15 


320 






922 79 


330 






983. 51 


340 






1047. 20 


350 






1112.31 


360 




6. 9629 
7. 1553 
7. 3421 
7. 5293 




1268. 71 
1339. 47 
1411.78 
1486. 17 




6. 7999 

6. 9939 
7. 1865 

7. 3798 




1179 27 


370 






1248. 29 


380 






1319. 13 


390 






1390. 97 


400 





7. 7141 





1562. 37 




7. 5727 




1466. 73 


410 




7. 8966 
7. 9748 
8. 0761 




1640. 42 
1674. 52 
1720. 26 




7. 7654 
7. 8482 
7. 9578 




1543. 41 


414.3 






1576. 98 


420 






1622 03 



















Liquid 


phase 








200 


a 6. 8 


6. 6360 


1076 


1064. 60 


4.83 


4. 8676 


-110 


-91.08 


210 


6.0 


6. 7178 


1144 


1131. 37 


5.16 


5. 1934 


-58 


-40. 77 


220 


6.95 


6. 7996 


1203 


1198. 96 


5.50 


5. 5078 


4 


12.75 


230 


7.0 


6. 8814 


1274 


1267. 36 


5.80 


5.8118 


62 


69.35 


240 


6.4 


6. S632 


1343 


1336. 58 


6.10 


6. 1064 


124 


128. 95 


250 


7.05 


7. 0450 


1405 


1406. 62 


6.37 


6. 3923 


187 


191. 44 


260 


7.0 


7. 1268 


1475 


1477. 48 


6.65 


6. 6702 


253 


256. 76 


270 


7.25 


7. 2086 


1544 


1549. 16 


6.92 


6. 9407 


324 


324. 83 


280 


7.35 


7. 2904 


1615 


1621. 66 


7.19 


7. 2044 


398 


395. 57 


290 


8.4 


7. 3722 


1688 


1694. 97 


7.45 


7. 4616 


470 


468. 90 


300 


7.7 


7. 4540 


1770 


1769. 10 


7.74 


7. 7129 


547 


544.77 


310 


6.75 


7. 5358 


1845 


1844. 05 


8.00 


7. 9586 


628 


623. 13 


320 


7.7 


7. 6176 


1911 


1919. 82 


8.22 


8. 1992 


719 


703. 93 


330 


7.6 


7. 6994 


2000 


1996. 40 


8.45 


8. 4348 


794 


787. 10 


340 


7.6 


7. 7812 


2072 


2073. 80 


8.69 


8. 6660 


883 


872. 62 


350 


7.55 


7. 8630 


2152 


2152. 02 


8.90 


8. 8926 


968 


960. 40 


360 


7.95 


7. 9448 


2225 


2231. 06 


9.12 


9. 1153 


1062 


1050. 44 


370 


8.5 


8. 0266 


2304 


2310. 92 


9.34 


9. 3341 


1154 


1142. 70 


380 




8. 1084 


2392 


2391. 60 


9.56 


9. 5492 


1241 


1237. 12 


390 




8. 1902 
8. 2720 
8. 3538 
8. 3890 




2473. 09 

2555. 40 
2638. 53 
2674. 52 




9. 7610 
9. 6693 
10. 1746 
10. 2619 




1333. 68 


400 






1432. 32 


410 






1533. 06 


414.3 






1576. 98 


420 




8. 4356 




2722. 48 




10. 3768 




1635. 80 











a This column refers to temperatures 5° higher than indicated; e.g., 6.8 is the 
C p value at 205 °K. 



value, and at the same time the correction was 
desired to have as little effect as possible on the 
linear thermodynamic functions for the solid below 
300 °K where a good fit of the data was already 
obtained. 
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Figure 1. Predicted specific heats, C PJ of an ideal CH 2 - 
chain crystal and liquid as a function of temperature. 

The closed circles represent values calculated from n-paraffin data, samples of 
which are listed in table 1, and the solid curves were calculated from the equations 
in table 3. 
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Figure 2. Predicted enthalpies (H-H ), of an ideal CH 2 - 
chain crystal and liquid as a function of temperature. 

The closed circles represent values calculated from w-paramn data, samples of 
which arc listed in table 1, and the solid curves were calculated from the equations 
in table 3. 
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Figure 3. Predicted entropies, S, of an ideal CH 2 -c/iam 
crystal and liquid as a function of temperature 

The closed circles represent values calculated from w-paraffin data, samples of 
which are listed in table 1, and the solid curves were calculated from the equations 
in table 3. 
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Figure 4. Predicted Gibbs free energies (F-H ) of an ideal 
CH.2-chain crystal and liquid as a function of temperature. 

The closed circles represent values calculated from ri-parafTm. data, samples of 
which are listed in table 1, and the solid curves were calculated from the equations 
in table 3. 

Fortunately, a correction of the type implied here 
by the r/ -paraffin data has the same form as the spe- 
cific heat contributions due to several independent 
CH 2 vibrations which become important above 
250 °K [8, 9]. These vibrations are well represented 
by a summation of Einstein functions of the form 
C v =Y^{Rd\IT 2 )[e e ^ T l(e^ IT -\) 2 \. For our use, we have 

i 

preserved the form of this function but have applied 
it to C ? rather than C v and have made the further 
approximations that we can use a single term to 
replace the summation of several terms, and that 
0>>T in the range of interest. That is, we shall 
assume a correction to C P of the form, 



C p =(nRd 2 IT 2 ) exp (-0/T), 



(1) 



where the bar indicates an "effective" quantity 
which is expected to lie somewhere between the maxi- 
mum and minimum quantities in the above expres- 
sion for C v . The effective number of independent 
CH 2 modes is represented by rl, 11 is the molar gas 
constant, d=hvjk is the characteristic temperature, 
h is Planck's constant, v is the effective frequency of 
vibration, k is Boltzmann's constant, and T is the 
absolute temperature. 

It is convenient to replace the two parameters 
n and in eq (1) by two different parameters, 
F s and H s , the corrections to the free energy and 
enthalpy at T M needed to give the proper melting 
point and heat of fusion. This replacement is ac- 
complished by integrating eq (1) to find the corre- 
sponding H and F and then solving the equations 
H S =H at T= T M and F=F at T= T M for 6 and n. 
The result is n = (F s /RT M ) exp (HJF 8 ) and 6=T M H J 
F s . If we assume T M =4U:S °K[4] then one finds 
that F, must equal —48.80 cal/mole (calculated 
from the uncorrected AF at 414.3 °K) and H s must 
equal (1162.10 — Ah f ) cal/mole where Ah f is the heat 
of fusion at T M =414.3 °K. We have here chosen a 
value of AA/=1000 cal/mole as reasonable, which 
choice gives 0=1376 and i = 1.64. 

The resulting analytical expressions for the thermo- 
dynamic funt ions of the CH 2 -chain crystal and 
liquid above 200 °K are shown in table 3, and the 
calculated values are compared to the experimental 
values in table 2. 



Table 3. Analytical expressions for the thermodynamic functions of the ideal OYLi-chain crystal and liquid above 200 ° K, derived 

from n-paraffin data 



Function 


Solid phase (T>200°K) 


Liquid phase (T>200 °K) 


C p 

(cal/mole °K) 


0. 0139 T+0. 916+22. 308 X 10* uj T* 


0. 00818 T+5. 00 


H-Ho 
cal/mole) 


0. 00695 T2+0. 916T-60+162. 10?* 


0.00409T2+5.00T-99 


S 

(cal/mole °K) 


0. O139T+0. 9161nT-3. 94+(0. 117789+162. 10/ T)u 


0. 00818 T+5. OOln T-23. 26 


F-Ho 

(cal/mole) 


-0. 00695 T2-0. 916 Tin T+4. 856 T-60-0. 117789 Tu 


-0. 00409 T2-5. 00 Tin T+28. 26T-99 




Liquid-solid difference {1 


^>200 °K) 


(cal/mole °K) 


-0. 00572 T+4. 84-22. 308X10%/ T* 




AH" 

(cal/mole) 


-0. 00186 T2+4. 84T-39-162. lOw 




AS 
(cal/mole °K) 


-0. 00572 T+4. 841n T+19. 32- (0. 117789+162. 10/ T)u 




AF 

(cal/mole) 


0. 00186 T2-4. 84 Tin T+23. 404 T-39+0. 117789 Tu 





w=exp[-3. 32172 (414. 3-T)/T] 
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6. Discussion 
6.1. Reliability of the Results 

It is important to have some idea of how accurately 
the results listed in table 2 represent the properties 
of an ideal CH 2 -cham compound. The specific heats 
below 150 °K and H, S, and F below 200 °K are 
probably accurate to ±1 percent judging from the 
consistency of the averaged values (examples in 
table 1) and the estimate of better than 1 percent 
accuracy for the original data. At 200 °K the effects 
of premelting can be seen in table 1 as a reduction 
in the first two specific heat differences and the first 
H, S, and F differences. One is increasingly un- 
certain above 200 °K about the correctness of the 
averaged data values, and thus the calculated values 
are to be preferred at the higher temperatures. All 
reasonable attempts to alter the results by starting 
with different initial linear and quadratic specific 
heat functions did not change the predicted melting 
point by more than 2 degrees or the predicted heat 
of fusion by more than 5 cal/mole. Nor did different 
assumed forms of the correction term make a 
significant difference in the predicted thermodynamic 
functions. The results were found to be quite in- 
sensitive to artifacts in the data or initial assump- 
tions. Tins insensitivity, while giving one confidence 
in the preducted values, greatly reduces the impor- 
tance which one can attribute to the values of 6 and 
n in eq (1). 

The most uncertain and sensitive parameter in 
the calculations is the heat of fusion, Ah f . The best 
current experimental estimate of AA f for polyethylene 
by the diluent method is 970 cal/mole [10]. A reason- 
able calculation shows that for the ideal CH 2 -chain 
crystal, Ah f can be expected to be higher than for 
100 A thick lamellar polyethylene by 2 to 3 percent. 
From our work, AA/— 975.9 appears to be a minimum 
value since an}^ lower value results in the solid spe- 
cific heat exceeding the liquid specific heat below the 
melting point. As Ah f is increased above its mini- 
mum value, the deviation between the calculated and 
experimental values in the 200 to 290 °K range in- 
creases, and the values of B and n in eq (1) decrease 
from 1580 °K and 2.70 for AA / =976 to 1370 °K and 
1.64 for A/^1000 cal/mole and 1164 °K and 0.98 
for Ah f =1025. If Ah f is as large as 1025 cal/mole, 
both the deviation from the experimental data and 
the values of 6 and n seem unreasonable. Thus it 
appears likely that the correct value of Ah f lies 
within ±2.5 percent of the 1000 cal/mole assumed in 
this analysis. 

The ±2.5 percent error in Ah f causes an uncer- 
tainty in H and S for the solid which decreases 
rapidly with temperature from about 1 percent at 
400 °K to a negligible error at 300 °K. The effect 
of the uncertainty in Ah f on the specific heat is 
illustrated in figure 5. 

For the liquid phase, the standard deviations of 
the data from the calculated curve have been given 
in the previous section and amount to roughly 10 
percent for C PJ 2 percent for F, and 0.3 percent for 
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Figure 5. The predicted specific heats, C p , of an ideal CH 2 - 
chain crystal and liquid (solid curves) for two assumed heats 
of fusion, compared to the extrapolated values given by Wun- 
derlich for completely crystalline polyethylene (dashed curve) 
and the polyethylene data of Dainton, Evans, Hoare, and 
Melia (closed circles) and Passaglia and Kevorkian (open 
squares) . 



II and S. Hence, it seems reasonable to expect the 
calculated H and S values to be accurate to within 
±2 percent. A comparison with polyethylene 
liquid data in the next section supports an estimate 
of ±2 percent reliability for the liquid C p values as 
well as the H and S values. 

6.2. Comparison With Experiment 

Figure 5 show r s a comparison of the ideal specific 
heat predicted here from n-paraffin data with that 
predicted by Wunderlich [1] from polyethylene data 
and with actual experimental values for crystalline 
Marlex 50 as reported recently by Passaglia and 
Kevorkian [9] and Dainton, Evans, Hoare, and 
Melia [11]. The specific heat values predicted here 
are appreciably lower than the polyethylene values 
up to 80 °Iv, slightly larger from 80 to 140 °K (this 
range varies somewhat with different data) and 
again smaller at higher temperatures. The smooth- 
ness (smaller variations in the slope) of the experi- 
mental C p curve relative to the predicted curve 
below room temperature indicates that the ideal 
vibrational frequency distribution is broadened and 
the characteristic temperature is lowered in actual 
specimens of polyethylene — a result attributable to 
imperfections. Similar behavior was noted during 
this study in C p data for a sample of n-C 32 H 66 and a 
less pure sample of 7i-C3sH 68 . Whereas the C 33 
specific heats were nearly linear with temperature 
from 90 to 180 °K, the C 32 specific heats w T ere 
markedly sigmoid al in the same temperature range 
with values below those of C 33 at high and low 
temperatures and above those of C 33 at tempera- 
tures between 100 and 150 °K. Liquid polyethyl- 
ene specific heat data [9], shown as squares in figure 
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.'). agree to within 1 percent with the predicted liquid 
C p curve. The enthalpies of the liquid reported by 
Passaglia and Kevorkian are smaller than predicted 
by less than 2 percent and the entropies of the liquid 
are larger than predicted by less than 2 percent. 
These agreements are considered satisfactory. 

Wunderlich 's predicted H and S values for the solid 
at 300 °K agree to within 2 percent with those pre- 
dicted from the ^-paraffins. The predicted specific 
heats differ somewhat as shown by the dashed curve 
in figure 5. Dainton, Evans, Hoare, and Melia [11] 
used an extrapolation technique to obtain a value for 
the ideal entropy of polyethylene at 25 °C of 5.5 
cal/mole °K in favorable agreement with 5.6 cal/ 
mole °K predicted here. 

One important result of this work concerns the 
thermodynamic crystallinity scale [i.e., the calcula- 
tion of the crystalline fraction of a sample from the 
equation x=[H (liquid) -H (sample) ]/[H liquid)-H 
(crystal)]. The values of H (liquid) -H (crystal) 
calculated in this study are significantly larger than 
those given by Wunderlich [1], but both the crystal- 
line and liquid enthalpies are smaller, probably 
because of a zero point enthalpy difference. Such 
a zero point enthalpy difference is expected if one 
assumes thai t he polyethylene extrapolation pertains 
to an ideally perfect chain folded lamellar crystal 
rather than the ideally perfect, large nonfolded 
crystal associated with the paraffin extrapolation. 
Whereas both types of perfect crystals would have 
the same zero point entropy (equal to zero by the 
third law), the zero point enthalpies would be differ- 
ent. It is interesting to note that the Marlex 50 
sample measured by Passaglia and Kevorkian [9] 
has a volume crystallinity of 80.0 percent and an 
enthalpy crystallinity of 93.8 percent as determined 
from the values of Wunderlich and Dole [12J, and 
using enthalpy values from this work, 97 percent or 
93.5 percent, depending upon whether one uses the 
solid at °K or the liquid above T M as the reference 
state. The zero point enthalpy correction which 
had to be added for the latter calculation was 40 
cal/mole of folds. The fact that thermodynamic 
crystallinities are generally higher than volume 
cryst affinities is consistent with the assumption that 
much of the imperfection in highly crystalline poly- 
ethylene samples is unoccupied volume (such as 
might be expected to exist between lamellar fold 
surfaces.) Volume measurements w r ould be sensitive 
to these defects whereas enthalpy measurements 
would not. 

6.3. Observations and Suggestions 

The thermodynamic properties represented in 
figures 1-4 exhibit several unusual features that 
result from the high specific heat of the solid at high 
temperatures (nearly equal to that of the liquid at 
the melting temperature). The difference (AH) 
between the enthalpies of the solid and liquid is 



quite constant in the melting region and the tendency 
of AH to become negative at low temperatures is 
not pronounced. Even though AS appears to tend 
toward zero between 100 and 150 °K, the necessity 
for polyethylene to form a glass to maintain a 
positive AH and AS is not as obvious as in most 
substances. It seems likely that the discontinuity 
in the AH versus temperature curve at the glass 
transition would be so small that it would be difficult 
to detect. In accord with the above, AF is well 
approximated by Ah f AT/T M over an unusually large 
range of temperatures. This simple linear approxi- 
mation to A/ is in error by only 1.5 percent for AjT=50 
degrees and by less than 5 percent for AT=100 
degrees. One would expect A/ to be much less 
linear than this for most substances [12]. From 
the above observations, it is not surprising that 
there is very little indication of a glass transition 
in the thermodynamic data for linear polyethylene. 

There is no indication in the data derived from 
the liquid paraffins above 200 °K of the existence of 
the glassy state in polyethylene, and yet the glass 
is generally considered to he the stable amorphous 
state in polyethylene below about 2^7 °K [1]. One 
may conclude that the glass transition temperature 
is chain length dependent and that the glass be- 
comes stable with respect to the supercooled liquid 
above 200° only for chain lengths greater than some 
minimum value. The only way one could use t Ik 1 
//-paraffins for predicting the properties of the 
polyethylene glassy state would be to have data, on 
the glassy state of the //-paraffins. 

As shown in table 1, the incremental increase per 
CH 2 in the n-paraffin thermodynamic properties is 
not as constant in the liquid state as in the solid. 
This result seems consistent with the following 
argument. In the solid, even for chains with less 
than 10 CH 2 units, the intermolecular forces are 
strong enough compared to the intramolecular 
forces to isolate the motions of a central CH 2 unit 
from those of the chain ends, in the liquid the 
intermolecular forces are weaker and less regular 
so that the CH 2 behavior is sensitive to chain end 
effects. The nonlinearity of the incremental increase 
per CH 2 in the liquid thermodynamic properties 
with reciprocal chain length possibly indicates the 
presence of "entire chain" vibrational modes in the 
liquid even for very long chain lengths. 7 

The departure of the predicted specific heat from 
that measured for polyethylene in the 10 to 80 
°K temperature range should be investigated further. 
Since all polyethylenes so far measured seem to 
show the same higher than predicted specific heats 
in this range, the anomaly would seem to be asso- 
ciated with some nonideal feature which all samples 
have in common — e.g., chain folded surfaces. It 
would be interesting to examine the changes in 



" In a private communication, ,T. P. McCullough lias suggested that the non- 
linearity in the liquid state is due to a contribution from short range order which 
decreases exponentially with temperature above the melting point. If this 
were the case, the liquid data would become more linear al higher temperatures. 
Higher temperature specific heat data can clarify this point. 
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specific heat of polyethylene with lamellar thickness 
to see if the discrepancy could be resolved. It 
might also be useful to attempt to analyze the data 
in this paper in the manner reported recently by 
Wunderlich [8]. It is desirable to have more n- 
paraffin specific heat data for longer chain lengths 
and higher temperatures in order to increase the 
accuracy and temperature range of the predicted 
higher temperature solid and liquid thermodynamic 
values. 
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